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Abstract 

We consider the cooperative decay of incoherently pumped atoms 
in a disordered medium, where light undergoes multiple scattering. It 
is shown that the cooperation number, which determines the duration 
and amplitude of super fluorescent impulses, is given by the number 
of atoms along a diffusive trajectory of the light propagating through 
the medium. We also consider the problem of reflection of a probe 
wave during cooperative emission. 

1 Introduction 

There is growing interest active photonic paints. These are media in which 
light undergoes multiple random scattering, resulting in a diffusive propa- 
gation of radiation, while interacting with atoms that can be pumped to 
obtain a positive population difference. The reflection and transmission of 
the electromagnetic waves through such a cavity has been extensively studied 
over the past decade. The speckle pattern resulting from scattering has an 
average enhancement in the direction opposite the direction of the incident 
radiation (a comprehensive review of other statistical properties of the 
speckle of reflected and transmitted waves is given in@). 

Feedback provided by scattering in such a random cavity can serve to set 
up laser oscillations || . The laser action in a powdered laser materials [|J || , 
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laser dye solutions with scattering nanoparticles ||, and dye-doped micro- 
droplets containing Intralipid as a scatterer |7| has recently been reported. 
These experiments concentrated mostly on temporal and spatial properties 
of emission. 

Recently, the proposed §| enhancement of the weak localization peak in 
backscattering from an amplifying photonic paint was observed ||. 

The relevant question concerning recent observations of generation of light 
in active photonic paints 0] ||is to what extent this phenomenon is remi- 
niscent of superfluorescence |I| |]TTJ| (i.e. , the cooperative decay of an inco- 
herently pumped system of dipole transitions, started by initial noise or an 
external electromagnetic field), which has usually been studied in systems 
without scattering. 

Here we consider the cooperative decay of incoherently pumped atoms in 
a random cavity, which is a slab of thickness L(L ^> I, where I is the mean 
free path of radiation). This geometry is often used in experiments. The time 
that light spends in this cavity is of order where D is diffusion constant. 
This time is to be compared with the energy exchange time between atoms 
and field. We show that if the latter is greater than then after some 
delay, the system will generate a superfluorescent pulse of hyperbolic secant 
form. 



The duration of the superfluorescent pulse is T ra dN^ []n| , where r rad is 
the time of radiative decay of a single atom and Nc is the cooperation num- 
ber, i.e., the number of atoms that take part in cooperative decay. We find 
that in disordered systems, this number is Nc ~ P^~t~ ( ^ is the wavelength 
of the radiation and p is the density of active atoms, such that pA 3 ^> 1), i.e., 
it is equal to the number of atoms in a tube with cross section A 2 and length 
of the order of j- , which is the length of a diffusive trajectory of radiation. 

The intensity of radiation of cooperating atoms at the maximum of the 
superfluorescent pulse is ~ — — JI0|. We show that the diffusive slab radiates 
at maximum as a system of independent groups of cooperating atoms, 
and at the peak of superfluorescent pulse, the intensity emitted by the slab 

AT 2 1/ 

is ~ — — x V is the volume of the slab. 

Trad 

The maximum cooperation number for given T rac i and density of active 
atoms is determined by the condition that the time of energy exchange be- 
tween atoms and the field equals the time that light spends in the cavity. 
From this condition, we find that the maximum cooperation number in the 
random cavity is NQ ax ~ \^/cpr ra ^ where c is the speed of radiation in the 
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slab. 

These results are valid in the case of weak dephasing processes and long 
relaxation of population difference. Below we take into account the effect of 
dephasing on superffuorescence. 

In the limit of the large escape time of radiation, atoms exchange 
energy with the field many times, so stimulated emission becomes important 
and the system exhibits oscillatory behavior. 

We also consider the reflection of the probe wave during decay of the 
pumped system. 



2 Basic equations. 



We model a random medium in the following way. The dielectric function 
e ("r*) of the medium, which contains active atoms, is a random function 
of position, such that (e(~r*)). Scattering of light is due to fluctuations 

of the dielectric function with white-noise like variance ( 8e ( 
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Se(7) 



We consider the case of a weakly disordered system / ^> A, with dimen- 
sions larger the than mean free path, so propagation of the field can be 
described as a diffusion process with diffusion constant D = c is the speed 
of light in the medium. 

The coupling between the polarization density | {e lujt P (~r*; t) + e _l£J *P* (V; £)}, 
averaged over scales smaller than A, the population difference density AN (~r*; t), 
and the field § {e iujt E (~r>; t) + e 
sical Maxwell-Bloch equations. 

emission noise is neglected, which is a good approximation for superffuores- 
cence [?]. P(T > ;t)and E(~r*;t) are slowly time- varying complex quantities, 
which we consider to be scalars; u is the atomic frequency. 



-iwt£j* ^-p>. ^ j can k e described by the clas- 
In this approach, amplified spontaneous 



First two Maxwell-Bloch equations have the form [13 



d 

dt +1 



P{?;t) 



(1) 



^-AN (r>; t) = ~ {P* ( r>; t) E (r>; t) - P (r>; t) E * (J>; t)} (2) 
at 2n 



Here 7 is the inverse dephasing time and \x is the electric dipole moment. 
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It is assumed that the population inversion relaxation time is longer than 
the delay time of the superfuorescent pulse. We also neglect inhomogeneous 
broadening. 

The quantities AN ( r ; t) and - ' are components of the local Bloch 
vector averaged over scales smaller than the wavelength of the radiation. The 
rate at which its length decreases, according to (1) and (2), is determined by 
T l . 

The field wave equation for the slow time-varying component E (7*; t) 
has the form 

dE(~r*:t) { c 2 u>e(j*)) ,_. . . 
—*-{-£;* ^\E(?;t) = 2™P(?;t) (3) 

Although E ("7*; t) and P (If*; t) vary slowly in time, they still contain 
spatial random phases, which result from random interference between waves 
coming to the point "r* via different diffusive trajectories. To get rid of 
these phase factors, it is convenient to consider the diffusion propagator 
D (7*; ti, £2), which determines the correlation function of the polarization 
density and field : 

(E ( r>; t x ) E* ( r>; £ 2 )> = 4ttA; 3 u; j D ( r>; t la t 2 ) (4) 

Correlation functions involving the polarization density can be obtained 
by using Eq.(l). 

To obtain the equation for the diffusion propagator, it is convenient to 
eliminate the polarization density from Eqs.(l) and (3). Then the usual 
diagram technique ]14| makes it possible to calculate the average of the 
product E (7*; ti) E* ("r*"; t 2 ) . 

Considering the evolution of the Bloch vector from time t = 0, at which 
a positive population difference is created, we obtain for the diffusion prop- 
agator 

1 h 

+ ^ [ dtexp{-j(ti -t)} AN(T y ,t)D(-r>;t,t 2 ) + (5) 

Q 
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1 h 

2 / rftexp{- 7 (t 2 -t)}A^(r,t) J D(r > ;t 1 ,t 2 ) 



Here To = \J 2np t\^\' 2 = y t^px^ * s ^ ne characteristic time of energy ex- 
change between the field and the atomic system [?], p is the density of active 

atoms, and r~^ d = is the radiative decay time of a single atom. 

The function / (V;^^) depends on initial conditions. Here we choose 
the initial condition such that (P ("f > , * = 0)) = and (p ("r > , i = 0) P* (7", t = 

p | /x | 2 5 ("r^ — r ). This initial condition corresponds to an initial incoherent 
state. In this case 

/ *i, *2) = P H 2 ex P {-7 (*i + *2)} (6) 

for times greater than the mean free time of radiation -. 

The equation for the mean population inversion density can be obtained 
by using Eqs. (1), (2) and (5): 



dAN (~r*,t) (2tt) 



dt hp\ 3 TQ 



J dtexp{- 7 (t-lf)}ANC?,lf){DC?;lf,t)+D(?;t,lf)} 



(7) 

For the population difference we choose AiV (T*, t = 0) = AiV > as the 
initial condition (AiV = p). 

The usual boundary conditions for the diffusion propagator are D ("r*; ti, t 2 ) = 
on an open surface and "r? V* D ("f > ; ti , t 2 ) = on a reflecting surface; "r? is 
normal to the reflecting surface. 

The diffusion approach is justified if the time of energy exchange between 
atoms and field is greater than mean the free time of radiation, r -. 



3 Cooperative decay in photonic paint 

Below we consider a slab of thickness L (L ^> I). Let z be the coordinate 
across the slab L y z >z 0. It is convenient to study the solution of Eq. (5) 
in the form 
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/L 00 
-J2^n(z)D n {t 1 ,t 2 ) (8) 
1 n=l 

where \l/ n (z) = sin ^ is an eigenfunction of the diffusion equation with 
boundary condition ^ n (z) = at the free boundaries z — 0, L. 

Let us consider the initial evolution of the diffusion propagator, when the 
population difference does not depend on time. For the coefficients in (8), 
we obtain from (5) 

^ d + 4r +uj n \ D n (t 1 ,t 2 ) = ^■p\/j,\ 2 Jdz^ n (z)exp{--f(t 1 + t 2 )} + 



dti dt 2 



AN 1 
Pr 2 o 

AN 1 
pr 2 



J dtexp{- 1 (t l -t)}D n (t,t 2 )i9) 


J dtexp{-j (t 2 -t)}D n (t u h) 



Here u n = ^jt- is an eigenvalue of the diffusion equation. Solving Eq. 
(9) via the Laplace transform with initial conditions D n (0, t 2 ) = D n (ti, 0) = 
(the field vanishes at t — 0), we obtain 



iMt,t)~exHU/R- 7 ) + (10) 



'u n \ 2 4A7V u T , 

The critical value of positive inversion density AN n , above which the 

growth rate of a particular diffusion mode Z n = \j — 7) + — ^ — 7 

becomes positive, is AN n = ^rprf [?]. More detailed calculations of (10) 
are given in the next section. 

To proceed further in solving Eqs. (5) and (7), we make two approxima- 
tions. 

1). Below we consider the case of fast escape of radiation from the system, 
where oj\ 3> 5^,7 (or, according to (10), tqU\ > 1 for weak dephasing), 
so we can neglect the time derivative in Eq. (5). In the language of su- 
perfluorescence, this situation corresponds to the case in which there is no 
energy exchange between the emitted field and atomic subsystem [?]. The 
field serves only to develop correlation between atoms. 
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2). We consider only the most unstable mode D 1 (ti,t 2 ). At t — 0, 
the off-diagonal elements of AN nm = j dzAN (z, t) \I/ n (z) \& m (z) are zero by 
definition, and interaction between modes is irrelevant for most of the time 
of decay. We therefore assume that the interaction of the first diffusion mode 
D\ {ti,t 2 ) with higher modes does not qualitatively change the description of 
cooperative decay. 

Under these assumptions the equation for the diffusion propagator has 
the form 



Di (t u t 2 ) 



71 '00 'i 



p|^| 2 exp{-7 (ti + t 2 )} + 




o 



and for the population difference 




(12) 



Introducing 



Di (h,t 2 ) 



exp{- 7 (t 1 +t 2 )}G(x (*i);x(*2)) 



(13) 



in (11), where 




o 



(14) 



we obtain 




(15) 



Equation (15) can be solved by a Laplace transform as 
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G{Xi;X2) = / tt-^— 16 

txuj\ J 2m 2m Z\Z2 — Z\ — Zi 

C—ioo 

The asymptotic form of (16) for xi — Xi = X > 1 is 

Gtex) ^2pH!exp^} 



7ru;i Jt^X 



(17) 



The equation for the population difference (x > 1) is 



> X = -7^— 2 G (x; X) exp {-2 7 t} (18) 



rf^ 8fc_ 3 

Taking into account only exponential factors, we obtain the solution of 
Eq. (18) : 



f 2SN(t -t) \ | 7 PT V 
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Here we introduce 5N = AN — lpT ° x ; AiV is the population difference 
at the beginning of exponential growth of radiative intensity, when deviation 
from the initial population difference is small (AiV = p). 

The delay time in (19) is to = ^t^" m {^"V'p^^x}- I* 1 deriving this 
expression we took into account the relation between the time of energy 
exchange between atoms and field and \p\ , which enters into the initial 
condition for polarization density. 

The radiative intensity is proportional to rfA J^ n , and is emitted as a 
hyperbolic-secant pulse. The result (19) coincides with that of the Markov 



theory of superfluorescence in a system without scattering |T0|, ?]. The dif- 
ference is in definition of the cooperation number. 

It follows from Eq. (19) that in the case of weak dephasing, the duration of 

2 

a superffuorescent pulse is = r ra dNQ , where T ra d is the time of radiative 
decay of a single atom and N c = 18 ^ L is the cooperation number, i.e., the 
number of atoms that take part in the cooperative decay. This is equal to the 
number of atoms in a tube of cross section A 2 with the length of the diffusive 
trajectory -p The maximum of the cooperation number is determined by the 
condition t co>i = 1, whereupon N™ ax = 2\y/Qcpr rad . Under this condition, 
atoms can exchange energy with the field only once, i.e., stimulated emission 
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can be neglected. We note that for a given density p, decay time r rad , and 
velocity, the maximum cooperation number in a disordered system is smaller, 
than in a pencil-shaped system without scattering |T5[ . 

The maximum emitted radiation is V-^ANn (t = to) {V is the volume 

of the slab). It can also be written Nx-^, where N= is the number of 

' T rad ' N c 

cooperative regions in the slab. The cooperative decay in a diffusive medium 
can therefore be interpreted as the independent cooperative decay of N= 
systems, each consisting of Nq atoms. 

Dephasing processes increase the duration of a pulse by the factor £—3 — 

AN— 7pT o" 1 

and decrease the peak intensity by the square of this factor. Note that this 
result coincides with that for a system without scattering [?]. 

If 1 ^> TqUx, atoms exchange energy with the field many times. In this case 



we expect spiking of intensity. The frequency of spiking can be estimated [13 



from Eq. (10) as y ^pf — ~ l) - To obtain this expression we insert a 

negative value of the population inversion [I3| into (10) (this situation will 
occur after the pumped atoms exchange energy with the field). 



4 Amplification in the backward direction 

Correlation between pumped atoms can also be due to the external field, 
which stimulates emission in the forward direction in a system without scat- 



tering |T3[. In a disordered system one might expect enhancement of emission 
in the backward direction. 

Here we consider the reflection of a weak probe plane wave with frequency 
during the development of superfluorescent emission. The amplitude of the 
probe is low, so the effect of the external field on emission can be neglected. 
We can also neglect interference between the external field with the emitted 
one, because the initial state of polarization is incoherent. Below we consider 
in detail the linear stage of decay when the inversion density is high enough 
to produce only the lowest diffusion mode instability, AiV = AiVi (1 + 5) , 
5< 1. This situation resembles the experimental setup of |J. 

It is convenient to calculate the albedo, which is the ratio between the 
intensities of the reflected and incident fields. The time-dependent albedo 
can be expressed as fTH 
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a{!f;t) 



oo 

C 



4tt/ 2 





J dzdz' exp | i — | y d"^ {1 + cos if 7^} D (z, z' , ~~f?; t, t) 



(20) 

Here if is the sum of the incident and outgoing wave vectors, and is 
the position in the plane. Diffusion propagator (20) obeys Eq. (5) with the 

substitution of 5 (~r* — r ) for / ("r*; t\, t 2 ). We also assume that the incident 
wave is close to the normal to the surface. 

The first term describes diffusion scattering, and the second term de- 
scribes the interference part, which is strongly peaked in the backward di- 
rection. The physical mechanism of the interference contribution is exhaus- 
tively discussed in the literature; see, for example Refs.. |], H and references 
therein. 

The diffusion propagator can be represented as 
D (V, 7\ t u t 2 ) = E *n 0) (z 1 ) exp {ilf ( 7? " 7) } D n (?; h, t 2 ) (21) 

n 

The Laplace transform of Eq. (5) for time-independent AiV > yields 

ioo+C , . 

' , . . n f dpidp 2 exp (pitt + p 2 t 2 ) 

" (?; " ,= -lc («) 2 ^{p 1+ p, + a,fe)-^fe + ^)} 

(22) 

where the real part of the integration contour passes to the left of all 
singularities, and Q n (q) = Dq 2 +u n is the eigenvalue of the diffusion equation 
for the slab geometry. 

At t\ = t 2 , integrating over the difference pi — p 2 in (22), we obtain for 
the first mode at t > u>i 1 

.3 ioo+C 

Dn / t t x = ( 2 7) 2 f dp_ exppt 

fafa + ty-Mo 2m Py/P ~ Zi (?) (p + ^g^^P ~ Z 1 (q) 

(23) 

Here we introduce Zi (g) = — 1^~) ' which is the growth rate of 

D n (q; t,t). 
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For moderate times < ^> we obtain from (23) 

D (a-tt) = 2{eXp{Zl{q)t} - 1) (24) 

This expression is valid for either sign of Zi, i.e., above as well as below 
threshold. 

Taking into account that ^1(2) = -y/^sin^ , we obtain the singular 
contribution to the albedo from the first mode : 

3/ ( exp {Zi (0) t} - 1 exp{Zi(g)t}- . , 
(«. *) = ^| 5 + > ( 25 ) 

Below threshold the albedo is saturated. The peak at large times has a 
laplacian form oc — Krr- At threshold and above there is narrowing of the 

peak with increasing time. Exactly at threshold the albedo is linear with 
time, and above threshold the albedo grows exponentially. 



5 Conclusions 

To summarize, superfluorescent emission of active photonic paint develops 
due to the cooperation of atoms along a diffusive trajectory through a system 
with cross-sectional dimensions of the order of a wavelength. The pulse there- 
fore becomes narrower with decreasing mean free path of radiation until the 
cooperation number reaches its maximum value. The maximum cooperation 
number does not depend on disorder. 

An external field enhances emission in backward direction. The peak 
sharpens in coherent backscattering during cooperative decay in a disordered 
system. 

We thank A.V. Gol,tsev for useful suggestions. This work was supported 
by the Russian Fund for Fundamental Research under Grant number 97-02- 
18078. 
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